Emergence of atomic-density waves in a trapped Luther-Emery fermion gas 
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We present a novel and comprehensive microscopic study of Luther-Emery-paired phases in a 
strongly interacting atomic Fermi gas inside a parabolic trap and a one-dimensional (ID) optical 
lattice. Our work is based on a lattice version of density-functional theory, which uses as reference 
system the ID homogeneous Hubbard model. We test our approach for repulsive interactions 
against Quantum Monte Carlo data and show that, for sufficiently strong attractions, an atomic- 
density wave (ADW) in the central portion of the trap breaks the discrete translational symmetry 
of the underlying lattice. We demonstrate that the emergence of an ADW has a dramatic impact on 
experimental observables such as the Fraunhofer diffraction pattern and the momentum distribution. 
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Strongly-interacting ID quantum liquids are nowadays 
available in a varietyof laboratory systems ranging from 
carbon nanotubes [l| to semiconductor nanowires 3, 
conducting molecules |j| , and gases in optical lattices 4| . 
Chiral Luttingcr liquids at fractional quantum Hall 
edges also provide a beautiful example of a ID quan- 
tum liquid. The effective low-energy description of all 
these II? systems is based on a harmonic theory of long- 
wavelength fluctuations |6j due to the interplay between 
topology and interactions. In the most interesting exper- 
imental situations translational invariance is broken by 
inhomogeneities of various types, such as Hall bar con- 
strictions in the case of quantum Hall edges or trapping 
for ultracold atomic gases pfl . These strong external per- 
turbations induce new length scales causing novel physi- 
cal behaviors relative to the corresponding unperturbed 
model system. 

A powerful theoretical tool to study the interplay be- 
tween interactions and inhomogeneous external fields of 
arbitrary shape is provided by density-functional theory 
(DFT) 110, which is based on the Hohenberg-Kohn the- 
orem and on the Kohn-Sham mapping onto an auxiliary 
noninteracting system. Many-body effects enter DFT via 
the exchange-correlation (xc) functional, which is often 
treated by the local-density approximation (LDA) 0, 0] . 
The essence of LDA is to locally approximate the xc en- 
ergy of the inhomogeneous system under study with that 
of an interacting homogeneous reference fluid, whose cor- 
relations are transferred by the LDA to the inhomoge- 
nous system. Most of the applications of DFT to in- 
homogeneous electronic systems use the interacting ho- 
mogeneous electron gas (HEG) as the underlying ref- 
erence fluid. The exact HEG xc energy is not known 
but can be calculated to a high degree of numerical pre- 
cision with the help of Quantum Monte Carlo (QMC) 
data. However, there are a few interesting examples in 



the literature in which either the reference system is 
not the HEG or the auxiliary system of the Kohn-Sham 
mapping is not an assembly of noninteracting particles. 
In particular, inhomogeneous ID fermionic systems ap- 
pear as an example in which it is appropriate to change 
the reference system to one that possesses ground-state 
(GS) Luttinger-liquid rather than Fermi-liquid-type cor- 
relations (see Ref. [ll| and references therein to earlier 
work) . 

In this Letter we present the first such DFT study of 
trapped ultracold Fermi gases. We focus our attention 
on a two-component Fermi gas in a II? optical lattice 
with cither repulsive or attractive inter-component inter- 
actions and in the presence of a static external potential. 
Creating these types of systems experimentally seems to 
be within the reach of present-day technology [12| . Stud- 
ies of this model have only been carried out in the case of 
repulsive interactions 0, IT3. Hfil ] whereas the main focus 
of this Letter is on attractive interactions. We use a lat- 
tice version of DFT within an LDA based on the II? ho- 
mogeneous Hubbard model (HHM) 0, which transfers 
Luttinger-Mott and Luther-Emery-type correlations 
to the inhomogeneous gas. A novel type of ground state 
is found at strong attractive coupling. 

We consider a Fermi gas with Nf particles hopping on a 
II? lattice with N s lattice sites. Each site is labeled by the 
discrete coordinate {zi =i,i E [l,N s ]}. We assume that 
the system can be described by a single-band Hubbard 
Hamiltonian |l8j . 

H = [ct(z i )c a (zj)+K.c.] 

+ U/] ^(zj)hi(zi) +yy f Mt,{zi)h(z i ) . (1) 

i i 

Here t^j — t > if i,j are nearest-neighbor sites and 
zero otherwise, a =|, \ is a pseudospin-1/2 degree of 
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freedom (hyperfine-state label), n CT (zj) = c\(zi)c a {zi) 
is the pseudospin-resolved site occupation operator nor- 
malized to the number of particles with pseudospin a, 
(Z)i^ff(^i)) = N a , and h(Zi) = J2a™<r( z i) is tne total 
site occupation with (X^"-( z i)) = -^f- Finally V cx t(z) is 
an external potential that we take later below as repre- 
senting a harmonic trap. 

If Vext{z) = the Hamiltonian in Eq. (JIJ reduces to 
a ID HHM that can be exactly solved using the Bcthc- 
Ansatz (BA) technique for both repulsive (U > 0) and 
attractive (U < 0) interactions 16] . The properties of 
the system are determined by three parameters: (i) the 
fillmg factor n = N{/N s , (ii) the pseudospin polarization 
s = (iVf — Ni) /(2N S ), and (iii) the dimensionless cou- 
pling constant u = U/t. The energy e(n,s;u) per site 
can be calculated by solving a set of nonlinear coupled 
integral equations. In agreement with the Lieb-Mattis 
theorem [l9j the GS has s = 0. In what follows the GS 
energy at s = will be denoted by s GS (n, u) — e(n, 0; u). 
For u > the ID HHM describes a Luttinger liquid |20j 
if n 7^ 1,2: the two metallic GS branches for n < 1 
and n > 1 are connected by particle-hole symmetry, 
£ GS (n,u) = (n — 1)U + e GS (2 — n,u). At half filling, 
i.e. for n = 1, the GS is a Mott insulator for every 
M 7^ 0, as signaled by a cusp in e GS and by a gap in 
the charge excitation spectrum, while for n — 2 the GS 
is a band insulator. Turning to the case u < 0, the 
ID HHM describes a Luther-Emery liquid |l7j ]. which 
exhibits a gap As(n) in the spin excitation spectrum 
and no charge gap. Due to an exact symmetry of the 
ID HHM, the GS energy for u < is given by the 
energy of a repulsive ID HHM at half filling and vari- 
able spin polarization s(n) = (1 — n)/2, e GS (n, — \u\) — 
— \U\ n/2 + e(l, s(n); \u\). The Bardeen-Cooper-Schricffcr 
(BCS) approximation has been shown \2l\ to provide a 
good description of the GS energy but a rather inaccu- 
rate and even qualitatively wrong description of the gap 
A s (ra). 

In the presence of a confining external potential H. 
cannot be diagonalized exactly. We calculate the GS 
properties of TL in Eq. JTJ by a lattice version of DFT, 
the so-called site-occupation- functional theory (SOFT). 
SOFT was introduced in a pioneering paper by Gunnars- 
son and Schonhammer [22] in order to study the so-called 
band-gap problem in the context of ab initio theories 
of fundamental energy gaps in semiconductors and in- 
sulators. Within SOFT the exact GS site occupation 
^gs(^z) = (GS|n(zj)|GS) can be obtained by solving self- 
consistently the Kohn-Sham Schrodinger equations 

y^[-tj,j + v KS (zi] [n GS (zi)])Sij](p a (zj) = e a ip a {zi) , (2) 

3 

with v KS (zi; [n GS (zi)}) = Un GS (zi)/2 + v KC (zi) + V eKt (zi), 



together with the closure 

^gs(2j) = T a \(p a (zi)\ 2 . (3) 

CK,OCC. 

Here the sum runs over the occupied orbitals and the 
degeneracy factors r Q satisfy the sum rule J2 a ^« = 
Nf. The first term in the effective Kohn-Sham poten- 
tial is the Hartree potential, while v xc (zi; [n as (zi)]) — 
5E xc [n(zi)]/6n(zi)\ GS is the xc potential defined by 
the functional derivative of the xc energy E KC [n(zi)] 
evaluated at the GS site occupation. The total 
GS energy of the inhomogeneous system is given 
by E GS [n GS (zi)} = J2 a ^a£ a ~ J2i v ^(zi)n GS (zi) - 
E. i C/^ s (^)/ 4 + ^xc[n GS (z 4 )]. 

Equations J2J and © provide a formally exact scheme 
to calculate n GS (zi) and E GS [n GS (zi)], but E KC and u xc 
need to be approximated. The LDA has been shown to 
provide an excellent description of the GS properties of 
a variety of inhomogeneous systems Q . In the following 
we employ a BA-based LDA (BALDA) functional [HI, 

< C ALDA (^; Mzi)]) = v% m (n,u)\ n ^ nGs{zi) . (4) 

Here the xc potential of the ID HHM is defined by 
v^° m (n,u) = d n [e GS (n,u) — k(n)] — Un/2, with k(n) — 
— At sin (n7r/2)/7r being the kinetic energy of the nonin- 
teracting gas. We should point out two appealing fea- 
tures of Eq. (0J from the formal DFT viewpoint. First 
of all, the xc potential of the reference ID HHM can 
be exactly calculated from BA integral equations [23j] . 
Secondly, the Mott cusp in £ GS (n, u) for repulsive inter- 
actions is responsible for a discontinuity A xc (f7) in v^° m 
at n — 1, A xc ({7) = U — 21im„_ 1 - d n e GS (n < l,u). As 
a consequence and contrary to the HEG-based LDA, the 
xc potential in Eq. naturally possesses a derivative 
discontinuity dij. 

We demonstrate the accuracy of the BALDA scheme 
for trapped ID Fermi gases in Eqs. ©-0 by showing 
in Fig. n some illustrative results for the GS site occupa- 
tion of harmonically trapped repulsive fermions against 
the state-of-the-art QMC data of Rigol et al. 13]. The 
agreement between our theory and the QMC results is 
clearly excellent [2^]. We now turn to describe how the 
density profile changes when U becomes negative. In 
Figs. [3 and [3] we report our theoretical predictions for 
the density profiles of attractive fermions. We clearly see 
from these figures that the consequences of Luther-Emery 
pairing in a confined Fermi gas are dramatic. The GS site 
occupation is a finite-size density wave, i.e. an oscillat- 
ing function of Zi that reflects the tendency of atoms with 
antiparallel pseudospins to form spin-singlet dimers. 

In Fig. |3 we illustrate the role of the external harmonic 
potential. For a very weak confinement the central part 
of the trap is occupied by a small-amplitude ADW, but 
in a sufficiently strong trap a region of doubly-occupied 
sites is formed at the center of the trap. In Fig. [3] we 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 





V 2 = 2x 10— 't 



6 « 



20 



40 



o 



BALDA 
QMC 

BALDA • 

QMC O 

BALDA ■ 

QMC □ 



Vj = ix 10" 3 i 



4 6 

J HLfiiik'bimm 

60 80 100 



20 3D 41) fiO 60 



- u - -4 










.............. 



0.8 - 
0.6 - 



JlIIIIII! 





--- -J -J V -f V 


a *i *i I) i 1 i 





FIG. 1: Ground-state site occupation n GS (zi) as a function of 
site Zi for a paramagnetic repulsive Fermi gas with u = +6 
and in the presence of a harmonic potential V c ^t{zi) = Vi{zi — 
N s /2) 2 , with Vi/t = 2 x 10" 3 ,4x 10 -3 , and 6 x 10" 3 . In this 
Figure and in all the others we have chosen Nf = 30 and 
N s = 100. The results of the BALDA self-consistent scheme 
are compared with the QMC data of Ref. liH . 




FIG. 3: Ground-state site occupation n GS (zi) as a function of 
site Zi for a paramagnetic attractive Fermi gas in the presence 
of a harmonic external potential with Vi/t — 4 x 10~ 3 and 
for u = —0.5, —2, —4, and —20. The thin solid lines are just 
a guide for the eye. 



TABLE I: Ground-state and pair-binding energies for a para- 
magnetic attractive Fermi gas in a harmonic external poten- 
tial for u — —0.5, —2, —4, and —20. Note that we are able to 
handle also large systems with our SOFT scheme. 
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FIG. 2: Ground-state site occupation n G s(zi) as a function 
of site Zi for a paramagnetic attractive Fermi gas with u = 
— 2 and in the presence of a harmonic potential, with Vijt 
increasing from 10 -4 to 10 _1 . The thin solid lines are just 
a guide for the eye. In the bulk region (40 < Zi < 60) and 
for V 2 /t=10- 4 , n GS {zi) can be well fitted to an ADW of the 
form n GS (zi) = no + ni cos (2irzi/\ ADV/ + (p) with no = 0.43, 
ni = 0.12, Aadw = 4.70 (to be compared with the prediction 
Aadw = 2/n = 6.67 for the homogeneous gas (l^l. and ip = 
2.17. 



show the crossover from the weak-coupling BCS regime 
to a strong-coupling regime where the spin-singlet dimers 
are tightly bound. On increasing the magnitude of the 
atom-atom attraction the amplitude of the Luther-Emery 
oscillations increases dramatically, giving rise to a large- 
amplitude ADW with period A ADW ks 2. The tendency 
of the system to form stable dimers is also confirmed 
numerically by a calculation of the pair-binding energy 
£ P = E GS (N { + 2) + E GS {Nf) - 2E GS (N l + 1), which turns 
out to be negative (see Table HJ. 

We proceed to illustrate the impact of the above ADWs 



on some measurable quantities. The key point is that 
the real-space Luther-Emery ADWs with their intrin- 
sic periodicity A ADW ^ 1 break the lattice symmetry 
and transfer structure to momentum space at wavenum- 
bers k w fc A Dw = ±27r/A ADW , away from the recipro- 
cal lattice vectors G v — 2nv, v — 0, ±1, ±2, .... As 
an illustrative example we show in Fig. ^ the elas- 
tic contribution to the light-scattering diffraction pat- 
tern [i^l, i.e. the Fraunhofer structure factor S c \(k) — 
exp (— ikzi)n GS (zi)\ , in the case of a relatively 
strong harmonic confinement [note that S e i(G„) = 1 and 
S e i(k + G v ) = S e }(k)]. In this case, bulk ADWs with 
wavenumber fc ADW = ±7r characterize the GS density pro- 
file for a sufficiently strong attraction. The emergence of 
a satellite peak at k = ±7r for u < —2 is evident in Fig.0] 

For the same parameters employed in Fig. 0] the mo- 
mentum distribution function of an attractive Fermi gas 
possesses a satellite peak at k = ±7r, which could be ob- 
served in a time-of-flight experiment. In addition to light 
scattering and time-of- flight-type experiments, trapping 
in a ring-shaped optical lattice |27l | can also unambigu- 
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FIG. 4: The Fraunhofer structure factor S e \(k) as a function 
of wavenumber k in the range —2n < k < 27r, for varying 
it in a paramagnetic Fermi gas inside a harmonic trap with 
V^/t = 4 x 10~ 3 . For clarity we have limited the z values 
to the interval [0,0.3]. The satellite peak at k = ±7r has 
magnitude 0.014 for u — —2 and 0.28 for u = —6. The inset 
shows a pictorial description of the unit-cell doubling, which 
drives the emergence of the satellites. 

ously reveal the Luther-Emery phases. The response of 
the system to an Aharonov-Bohm flux [2^ can be used 
to characterize microscopically the nature of the ground 
state. In fact, if the GS has a spin gap as in the Luther- 
Emery phase, it has been shown [29j that for a large 
system the period of the GS energy as a function of the 
Aharonov-Bohm flux is ir, exactly half the expected one. 
Thus ultracold Fermi gases with attractive interactions 
could be a highly tunable system in which it is possible 
to observe direct evidence of Luther-Emery pairing. 

In conclusion we remark that our DFT scheme is also 
applicable to a number of other experimentally interest- 
ing systems such as non-paramagnetic gases, gases under 
the effect of time-dependent external potentials, and ID 
quantum electron liquids. 
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